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1 Introduction and Motivation 


Alternative versions of this topic might be, “do secondary particles retain some ‘mem¬ 
ory’ of the direction of the primary parent particles?” Of course they do in the beam 
fragmentation (large Feynman-a:) regions as studied in Fixed Target and ISR experi¬ 
ments. But what about for the central regions near rapidity [T] y = 0 for hadron col¬ 
liders? What mechanisms are at play for the forward-backward asymmetries, denoted 
Apb"^- What occurs in the transition between the central and beam fragmentation 
regions? 

The dehnition of Forward and Backward is usually pretty intuitive. Forward 
particles are those where the produced quarks (anti-quarks) or charge follows the 
direction of the beam particles with similar quarks, anti-quarks, or charge. The 
proton direction {y > 0) is forward for t-quarks (and jets), 6-quarks (and jets or decay 
particles), or baryons (A° and A°) for these studies), while the antiproton direction 
{y < 0) is forward for t-quarks (and jets), 6-quarks (and jets or decay particles), or 
antibaryons. Forward B~ mesons {hu) follow the proton direction, while forward 
(bu) mesons follow the antiproton direction. At the pp LHC, all baryons would be 
considered “forward” and antibaryons “backward” regardless of their rapidity. 

The forward-background asymmetries for the production of tt [2] m and W± B! 
have been extensively studied for pp collisions at the Fermilab Tevatron. Unlike the 
symmetric pp initial state at the LHC, the pp state has the beam valence quarks 
traveling in the initial p direction (y > 0) and the beam valence antiquarks traveling 
in the initial p direction (y < 0). Is this forward-backward asymmetry in the direction 
of motion of the initial state quark charges retained for the produced particles? 



Figure 1: Graphs for pp production of 66 pairs. For pp tt pairs, the qq annihilation 
graphs dominate. 

Figureshows the hrst order pp production graphs for the central, small | r/ |, re¬ 
gion for gluon fusion (dominant for 66) and qq annihilation (dominant for tt). Neither 
of these hrst order graphs produce a forward-backward asymmetry. However, inter¬ 
ference with the loop and initial and hnal state gluon radiation graphs can produce a 
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significant forward-backward asymmetry. Recent measurements indicate that ApB{tt) 
~ -h 13%, significantly different from zero. These experiments have motivated the 
calculation of the Next-to-Leading Order (NLO) and Next-to-Next-to-Leading Order 
(NNLO) processes to attain agreement with the observations. 

For bb production in the central region, the gluon fusion graph is dominant and 
the higher order interference terms are small, leading to an expectation that Apsibb) 
would be small. 

The Fermilab Tevatron is a good place to study ApB since the symmetry of the 
initial pp state has net charge C = 0 and net baryon number B = 0. The production 
cross sections for particles and antiparticles are expected to be related as 

a{particle, y) = a{partide, —y) 

compared to the LHC where the initial pp state has C = +2, B = +2, and 

a{particle, y) = a{partide, —y) 

and where a{partide,y) is not necessarily equal to a{partide,y). Thus the pp 
Tevatron allows complementary set of forward-backward comparisons with a comple¬ 
mentary set of expected symmetries for the produced particles. 
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The DO Detector 



The DO experiment, illustrated in Figurej^ is, at least to first order, forward backward 
symmetric with respect to rapidity ± y or pseudorapidity [T] ± r^. 



Figure 2: The DO Experiment for Run II of the Fermilab Tevatron [5]. 
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Since the charged particle tracking is based on scintillating hbers and silicon strips 
whose performance is independent of the magnetic field direction, DO regularly re¬ 
verses the polarities of the tracking solenoid and muon toroid magnets which cancels 
many of the detector-related or detector-induced symmetries in the raw data. 

The DO results presented here are based on the full Tevatron Run II data set of 
10.4 fb~^ integrated luminosity for pp collisions at ^/s = 1.96 TeV. 


3 Analysis for Aps 

Prior Tevatron and LHC studies of the asymmetries for the ti and bb systems used jets 
tagged by the charge of the decay leptons, and used as their metric the difference in 
rapidities Ay = Pq-jet — Vq-jet- In this analysis, DO studies asymmetries for the single 
particle A°/A^, and A°/A^ inclusive production, without requiring the detection 
of both members of the B~ and , or A° and A°, or A° and A° pair in the same 
event. 

The Forward-Backward Asymmetry for the production of a particle at rapidity y 
is defined as 


. . X _ tTF(|/) - asiy) _ DIFF 

^ <JF{y) + (JB{y) SUM 


where (Tpip) Asiy)) is the differential cross section daF{y)/dy (dasiy)/dy)) for 
the Forward (Backward) production process. The asymmetry, as usually defined, is 
the difference divided by the sum for processes being compared. 

We measure the Afb asymmetry simultaneously for the particle and antiparticle 
and y < 0 and y > 0 hemispheres and define 


Afb(I y\) = 


apiparticle, y > 0) + apipartide, y < 0) — aB{partide, y < 0) — aB{partide, y > 0) 
(Tpipartide^y > 0) -|- <TF{partide,y < 0) -|- <TB{partide,y < 0) -|- aB{partide,y > 0) 


To assure no mixing of the forward and backward hemispheres due to reconstruc¬ 
tion errors, particles with | y | < 0.1 or | r] | < 0.1 are excluded from these sums. This 
assures that 99.9% of the have the same sign of y or p of its parent quark. For 
reference, the distribution of {pb-quark — VB-meson) has an r.m.s. width of 0.11 unit. 

The range of these Afb measurements is 0.1 < | y | < 2 units of rapidity. 


4 Afb{B^) |®i 

Charged B-mesons were reconstructed as B^ -d J/%l)K^ with subsequent decay 
J/V' —)■ //■'■/«“. Since DO does not have identification of hadronic particles, the charged 
decay product was assigned the mass of the kaon. The muons were each required to 
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have transverse momentum prifJ-) > 1-5 GeV and Pt{K) > 0.7 GeV. A Boosted 
Decision Tree (BDT) selection was used to help separate signal from background. 

An unbinned likelihood distribution was used to fit the invariant mass distribu¬ 
tions for the SUM = Forward + Backward and DIFF = Forward — Backward 
candidates including four contributions: 1 ) An exponential combinatorial background 
distribution; 2 ) A double Gaussian for the B^ J/%l)K^ state of interest; 3) A dou¬ 
ble Gaussian for incorrectly identified B^ —?• J/ipTr"^; and 4) An arctangent threshold 
function for partially reconstructed B^ —?• J/xph^X. For illustrative purposes, the 
Sum and Differences are plotted as binned distributions in Figure]^ There are 89,000 
B^ J/xpK^ candidates in the SUM distribution. A^b is the number of events in 
red DIFF distribution, divided by the number of events in the red SUM distribution 
for ^ J/xpK^. 


SUM = Forward + Backward DIFF = Forward - Backward 



Figure 3: Sum (Forward -|- Backward) and Difference (Forward — Backward) Mass 
distributions for B^. 

Integrated over 0.1 < | 77 | < 2, we find 

Afb{B^) = [-0.24 ± 0.42(stot) ± 0A9{syst)]% 

limited by the statistical uncertainty and consistent with no significant FB asym¬ 
metry. The ± 0.19% systematic uncertainty is dominated by variations due to alter¬ 
native boosted decision trees and cuts. 

An MC@NLO [7] plus HERWIG [ 8 ] simulation predicts — [—2.31± 

0.34(stat)±0.44(sr/st)] % for the same 0.1 < | 77 | < 2 range, which differs significantly 
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from zero asymmetry and is systematically higher than the DO measurement for all 
values of I r] |. An improved analytic calculation by C. Murphy of Pisa [H] gives 
= [0.021 ± 0.008]% for the same r] range, while reproducing the much larger 
observed ApB{tt). These measurements, simulations, and calculations are compared 
in Figure]^ as a function of | |. 


£ 


4 

3 

2 

1 

0 


DO data 
SM (this woik) 

1 MC«NLO+Hen.^g /WC@A/Z.O 


'-'>-1- i> -1 

J_I_ 1 I _I_I_I_I_l_l_I_I_I_I_ 1 _I_I_I_I_I_L 


0.5 1.0 1.5 2.0 


m)\ 


Figure 4: Comparison of measured ApBiB"^) with predictions of MC@NLO [Zj sim¬ 
ulation (band) and recent calculation by C. Murphy [H] (red line). 
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The decays A° —>■ J/'ipA^ with subsequent decay J/ip ^ yAiJi~ and A° pn~ (and 
the corresponding antiparticles) required a transverse separation of the A decay and 
primary interaction vertices of 0.5-25 cm, and a significance of the transverse sepa¬ 
ration of the A° decay vertex from the primary vertex of L^ylcri^y > 3. The average 
transverse momentum is < pt(A°) > = 9.9 GeV. Since DO does not have hadron 
identification, the higher momentum particle for the A^ or A^ candidate was assigned 
the mass of the proton. 

In Figure]^ the data were fitted with a binned maximum-likelihood method with a 
Gaussian signal plus a second-order Chebyshev polynomial background distribution, 
resulting in a total of 842 ± 48 A° -|- \ events. 

The asymmetry integrated over 0.1 < | r/ | < 2 for < Pt(%6 ) > = 9.9 GeV is 
measured to be Afb{^) = 0.04 ± 0.07{stat) ± 0.02{syst). 
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Figure 5: Mass distributions for the range 0.5 <| y |< 1.0 for Forward and Backward 
^ J/ijjK and ^ J/ipK candidates. 



Figure 6: Forward-Backward Asymmetry for A^ and A^, compared with prediction of 
the Rosner String Drag m and the Heavy Quark Recombination ra Models. 

As shown in Figure]^ ^ function of | y | is consistent with zero but 

does show a tendency to increase with increasing \ y \. This is consistent with the 
tendency of the String Drag model by Rosner m in Figure where the di-quark 
in the proton fragment tends to predominantly pick up the produced 6-quark giving 
an average increase in longitudinal momentum < Apz > = -|- 1.4 GeV for the 
compared to — 1.4 GeV for the A;,. In Figure the Heavy Quark Recombination 
Model of Lai and Leibovich im increases very little over this | y \ range and 
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is consistent with — 0 with the DO measurement. 




Figure 7: (left) The Rosner String Drag Model [H] and (center and right) the Heavy 
Quark Recombination Model of Lai and Leibovich [12] . 


Both CMS and LHCb observe the ratio of = a{Kb)/<^pp{Rh) and plot this 

ratio as a function of the rapidity loss relative to the proton beam yheam — Ukb- Using 
the expected pp charge symmetry a{Ab,y) = a{Ab, —y), it can be shown that for pp 
reactions that 


nLHC _ n _ _ 
^VV 


app^Backward) 

app{Forward) 


where 


Rpp{,Rbi 12 / 1 ) 


1 — ApsiRb, I y I) 
1 + ApsiAb, I y I) 



Figure 8: The DO measured values of Rpp = a{B)/a{F) for A;, compared to Rpp = 
a{Ab)/a{Ab) for CMS [I3| and LHCb [Tl|. 
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Thus the DO Ratio R = asiRb+^bj/o'Fi-^b+^b) can be directly plotted in Figure]^ 
along with the CMS [12] and LHCb [Hj data for R = a(Ab)/a{Ab) showing similar 
behavior as a function of the rapidity loss Ay = ybeam — UAt, ■ 

6 Afb{-^^s) 

The A° particles were collected under three separate trigger conditions; 1) prescaled 
beam crossing or minimum bias triggers (these are called the inclusive A data set and 
consists of 2.3 Million events); 2) events with a A or A along with a J/'ip which decays 
to ijA iJi~ (0.7 Million events); and 3) and events with a A or A along with either a 
single yA or /r“, regardless of the charge and particle/particle correlation between 
the A the A or A (53 Million events). The DO primary two-muon trigger resulted 
in the very high statistics for the J/^pA sample. We have not fully understood the 
correlations between the charge of the single pA and the A or the A, so we cannot yet 
describe the full physics import of the single yA events. 

The plain, old time A° was reconstructed via the decay A° —?• ptt~ and A^ —?• p7T~^ 
as shown in Figure]^ Again, since DO has no particle identification for baryons, the 
particle with the higher momentum in the decay is assigned the proton mass. The 
A candidates were required to have 2 < Pt{A) < 25 GeV. 
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Figure 9: (left) An example of a representative Mass spectrum for A ^ pTr~^ and 
A —7- pTv~ candidates, (right) Measured Afb as a function of | r/ | for inclusive A and 


A. 


DO has not submitted the A and A analyses for publication, so these Afb{-A) 
results must be considered preliminary. 
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Figure 10: Apsi^) as a function of | y | for inclusive A, J/'ipA, and /rA. 

For the inclnsive A and A sample, in Figure]^ Apb{\ y |) increases with increasing 
I y I and is statistically significantly greater than zero forl<|y|. Apb{\ y \) 
for the J/ipA and for the pA data sample also has similar behavior in Figure 
Similarly, when plotted as a function of the rapidity loss relative to the beam proton 
or antiproton Ay = ybeam — Ua, the comparison of A° and A^ prodnction at the pp 
Tevatron and for the pp LHC [161 HZl [E] ; for the pp ISR [HI |20] , and also for the 
p — Be and p — Pb Fixed Target data 

^ ^ <^Pp{A,y) ^ l-AfB{A,y) 
app(A,y) l + A%{A,y) 

exhibit very similar behavior in Fignre[^ which wonld hint at a common nniversal 
function linking the central production region near y = 0 which is dominated by hard 
production of quark-antiquark and particle-antiparticle pairs and the forward beam 
fragmentation regions at high Feynman-a: = PpartidelPbeam- 
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preliminary DO, 10.4 fb 



Figure 11: Comparison of the DO measured Rpp — a{B)/a{F) and Rpp = a{A)/a{A) 
for LHC and ISR experiments and RpA for Fixed Target experiments as a function of 
rapidity loss Ay = ybeam — Va- ATLAS [IE], STAR [U], LHCb [TH], R-607 [IHI, R-603 
PU] , and Fermilab E-8 [H]. 


7 Conclusions and Things to Remember about the 
DO Studies oi Aps 

Apsitt) is approximately +13% over the same rapidity range. The early difference 
between theory and observation motivated calculations of the Standard Model to 
Next-to-Next-to-Leading Oder (NNLO) to produce agreement with the data. 

Afb{BA consistent with zero. The MC@NLO simulation predicts an asymme¬ 
try of approximately +2%, while the recent calculation of Christopher Murphy gives 
approximately zero asymmetry. 

Ai7’B(A°) is also consistent with zero. More statistics are required to differentiate 
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between the rapidity dependence of the String Drag or the Heavy Quark Recombina¬ 
tion models. The DO measurement of R(A°) is consistent with the R{Ay) dependence 
for CMS and LHCb as a function of the rapidity loss Ay = ybeam — Ua- 

exhibits a statistically significant increase with increasing \ y \. As 
the rapidity loss Ay decreases, R(A°) seems to exhibit a transition from the central 
rapidity region at colliders to the beam fragmentation region previously studied at 
Fixed Target and ISR experiments. This behavior suggests that there might be a 
universal curve as a function of the rapidity loss Ay. 

DO continues to produce results complementary to those of the LHC, especially 
those unique to pp collisions. Stay tuned for more results! 
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